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ABSTRACT 

This  study  was  an  evaluation  of  a  radiational  scheme 
for  inclusion  in  a  numerical  weather  prediction  model. 
This  scheme  uses  empirical  formulations  for  atmospheric 
absorptivities ,  scattering  and  reflectivity  and  cloud 
and  earth  surface  reflectivity  to  compute  the  solar  in- 
solation received  and  absorbed  by  the  earth  and  key  atmo- 
spheric layers.   Longwave  cooling  by  the  earth  is  also 
treated  with  empirically  derived  emissivities  for  water- 
vapor  and  carbon  dioxide.   This  radiation  model  includes 
a  two-layer  cloud  parameterization  scheme  and  calculates 
the  effect  of  an  additional  cloud  layer  on  previous 
studies . 

In  testing  the  validity  of  this  model,  use  was  made 
of  Fleet  Numerical  Weather  Central  data  in  the  form  of 
vertical  soundings  on  a  spring  day  over  the  ocean.   This 
data  was  used  to  calculate  the  radiation  balance  at  the 
tropopause,  the  earth's  surface  and  the  key  layers  of 
the  atmosphere. 
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I.    INTRODUCTION 

This  study  was  an  attempt  to  produce  an  empirical 
radiative  parameterization  for  inclusion  in  the  Fleet 
Numerical  Weather  Central  (FNWC)  primitive  equation  pre- 
diction model. 

The  radiative  parameters  described  here  include  two 
cloud  layers  of  fractional  amounts  CL(1)  and  CL(2)  at 
each  gridpoint  instead  of  a  single  cloud  layer  as  pre- 
viously used  (Plante,  1973).   The  possibility  of  an  addi- 
tional cloud  layer  tends  to  add  substantial  complexity  to 
the  radiation  physics,  but  in  general  this  paper  repre- 
sents an  extension  of  that  described  previously  by  Martin 
(1972) .   Some  of  the  details  of  this  radiation  model  are 
also  based  upon  an  unpublished  manuscript  (Martin,  1974), 
which  in  turn  has  some  similarities  to  the  multi-cloud 
models  of  Arakawa  (1972)  and  Rodgers  (1967). 

The  vertical  scale  used  in  this  study  was  based  on 
the  five  sigma  (a  =  p/tt)  levels  of  the  FNWC  primitive 
equation  model.   The  levels  in  this  study  are  called  k- 
levels  (k  =  1,  „..,  10).   The  k-levels  are  used  for  the 
purpose  of  data  input  as  shown  in  Figure  1.   Cloud  cover 
is  dealt  with  by  means  of  a  two-layer  parameterization, 
with  one  high-level  cloud  layer  from  k  =  4  to  k  =  6  and 
a  low-level  cloud  deck  from  k  =  8  to  k  =  9,  in  amounts 
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Figure  1.   Five-layer  radiative  sounding  used  in  this 
study.   Levels  are  identified  by  their  values  on  the  k  - 
scale,  while  layers  are  identified  by  their  level  boundary 
indices  (e.g.,   (8,  10)  in  parentheses.   Pressure-scaled 
water  vapor  and  CO   mass  increments  Au  and  Ac,  respective- 
ly, as  well  as  mixing  ratio  (q)  are  introduced  at  the  odd 
levels  while  temperature  is  introduced  at  all  levels. 
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CL(1)  and  CL(2),  respectively.   In  the  computations  for 
pressure-scaled  water  vapor  and  C02  absorber  masses  a 
realistic  method  was  used  to  extrapolate  water  vapor  mix- 
ing ratio  values  above  the  level  of  k  =  4,  which  is  nor- 
mally the  highest  moisture  level  available  from  radiosonde 
data.   The  system  for  computation  of  the  net  flux  of  ter- 
restrial radiation  makes  use  of  new  empirical  absorptivi- 
ties  of  water  vapor  and  CO   (Sasamori,  1968),  which  are 
incorporated,  in  the  Martin  (1974)  manuscript. 

To  evaluate  the  radiative  heat  transfers  of  this  study, 
three  meridians  derived  from  FNWC  analyses  were  used. 
These  meridians  are  all  from  the  Pacific  Ocean  area  and 
most  of  the  radiative  soundings  were  in  the  Northern  Hemi- 
sphere . 
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II.   DATA  PREPARATION 

The  data  used  in  this  study  was  in  the  form  of  three 
meridians  of  25  vertical  soundings  over  the  Pacific  Ocean. 
These  meridians  were  located  at  125°W,  170°W  and  145°E  as 
shown  in  Figure  2. 

The  soundings  on  these  data  lines  were  taken  from  the 
original  FNWC  analyses  of  temperature  at  the  surface  and 
at  nine  other  standard  levels  up  to  100  millibars  for 
0000GMT,  25  April  1973.   The  six  moisture  values  were  given 
by  the  vapor  pressure  at  sea  level,  and  by  dew-point  de- 
pression for  levels  from  925  mb  to  400  mb .   Table  1(a) 
gives  as  an  example  the  original  sounding  at  gridpoint 
(1,  1)  in  the  FNWC  map.   The  data-processing  methods 
applied  to  the  temperature  and  mixing  ratio  data  for  the 
determination  of  the  water  vapor  and  CO-  absorber  masses, 
and  the  cloud  amounts  within  required  k-level  boundaries 
(Figure  1)  are  detailed  in  the  following  subsections.   It 
should  be  noted  that  all  soundings  start  at  sea  level, 
and  as  a  reasonable  approximation  may  be  considered  to 
have  surface  pressure  TV  =  1000  mb.   Hence  the  k-levels  of 
Figure  1  became  1000,  900,  800,  700,  600,  500,  400,  300, 
200,  and  100  mb,  respectively. 

A.    TEMPERATURE  TERMS 

The  temperatures  were  listed  at  each  mandatory  level 
between  1000,  ...,  100  mb.   The  radiative  sounding 
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Fiqure  2.   FNWC  polar  stereographic  grid  and  meridians 
9    ,  selected  for  study.   The  longitudes  A  and  hour 
angles  of  the  sun  relative  to  local  noon  are  shown  for 
each  meridian . 
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TABLE  1(a).   Form  of  original  FNWC  grid  soundings.   The 

humidity  parameter  between  925 400  mb  is  dewpoint 

depression . 


Pressure  (mb) 

1000 
925 
850 
700 
500 
400 
300 
250 
200 
150 
100 


T  (°C) 

28.5 
20.0 
16.2 
7.1 
-9 

-21 

-36 

-45 

-55 

-66 


6 
4 
4 
3 
5 
9 
-80.3 


Humidity  Parameters 

27.7  mb 
3.2°C 
2.9°C 
5.0°C 
4.7°C 
2.0°C 


TABLE  Kb).   Form  of  the  corresponding  radiative  sounding 
with  temperature  and  mixing  ratio  listed  at  k-levels. 


Pressure  (mb)  T  (°K) 

1000  299.0 

900  291.7 

800  286.6 

700  280.3 

600  272.6 

500  263.6 

400  251.8 

300  236.8 

200  217.7 

100  192.9 

0  192c9 


Mixing  Ratio 


18 
16 
12 

7 
4 


35 
50 
97 
84 
75 
2.97 
2.21 


g(Jcg) 


-l 


73 

23 

,04 

,00 
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temperatures  at  certain  k-levels  were  obtained  by  inter- 
polation.  The  temperature  T    was  set  equal  to  the  listed 
air  temperature.   At  the  top  of  the  atmosphere,  the  tem- 
perature was  assumed  to  be  isothermal  from  100  mb  to  p=0, 
with  the  mixing  ratio  also  approaching  zero  at  p=0. 

The  interpolation  scheme  used  was  the  three-point 
Lagrangian  formula  of  the  form: 


Tk  =  T0((Pk-Pl)(Pk-P2))/((P0-Pl)(P0-P2)) 

+  T1((Pk-P0)(Pk-P2))/((P1-P0)(P1-^2)J  (2-1) 

+  T  ((P  -P  )(P  -P  ))/((P  -P  )(P  -P  )) 
2V   k  0'   k  k   l///vv  2   0    2   1 

Here  P„  lies  above  P, ,  where  k  is  the  pressure  level  de- 
0  k'  c 

sired  for  T  -interpolation,  and  P  ,P   lie  below  level  k. 
For  example,  with  k  =  9,  P   and  P   were  taken  as  925  mb 
and  1000  mb ,  respectively.   The  parameters  T  ,  T   and  T 
are  temperatures  at  the  same  levels  as  P  ,  P  ,  and  P  , 
respectively. 


B.   MOISTURE  TREATMENT  IN  THE  SOUNDING 

The  moisture  parameters  from  the  original  sounding 

(Table  1(a))  were  converted  into  the  mixing  ratios  at  each 

of  the  original  sounding  levels.   At  k  =  10  the  surface 

vapor  pressure,  e  .  ,  was  used  to  calculate  the  mixing 

air 


ratio  by  use  of  Equation  (2-2) 


q „  =  .62197  e  .  /1000  (2-2) 

^1000  air 


where  e  .   is  in  mb. 
air 
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In  calculating  the  mixing  ratios  at  the  remaining 
original  data  levels  Equation  (2-3)  was  used  (after  p.  63 
Fleagle  and  Businger,  1963): 


q(p)  = 


(A-  £) 
.62197  e_ T^_ 


exp 


-B(T-T  ) 

[ — ] 

T  T_ 


1000      (2-3) 


Here  q(p)  is  the  mixing  ratio  at  level  p. 

A  =  21.656 

B  =  5418. 0°K 

T  =  temperature  at  level  p  (°K) 

T   =  dew  point  at  level  p  (°K) 
D 

T-T   =  dew  point  depression  (°K) 
The  q-values  were  then  interpolated  to  the  k-levels  using 
an  equation  similar  to  Equation  (2-1).   The  results  of 
this  interpolation  are  shown  in  an  example,  Table  1(b). 

Accurate  humidity  information  is  not  generally  avail- 
able at  or  above  300  mb .   To  find  the  q-values  at  k  =  3 , 
2  and  1,  an  extrapolation  formulation  was  used 


q(p)  .  (-£—)  X 
q„      ^500' 


(2-4) 


This  formula  is  similar  to  that  developed  by  Smith  (1966)  . 
However,  instead  of  using  the  surface  values  of  pressure 
and  moisture  in  Equation ( 2-4 )  use  was  made  of  the  pressure 
and  moisture  values  at  the  500  mb  level.   Equation  (2-4) 
was  solved  for  the  parameter  X    using  the  "best-fit"  between 
the  variables  y  =    log  |-   and  x  =    log  j^    .       This  was 
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done  by  using  the  six  values  of  q  of  each  original  sound- 
ings corresponding  to  p  =  1000,  925,  850,  700,  500,  and 
400  of  Table  1(a).   The  "best-fit"  values  of  X  were  then 
used  in  Equation  (2-4)  to  extrapolate  the  q-values  to 
p  =  300,  200  and  100  mb . 

The  usual  formula  for  a  best-fit  parameter  X  relating 
y   and  x.  samples  is  given  in  standard  texts  on  statistics 

(e.g.,  Crow  et  al,  1960) 

i=6 

E   *iXi 
i  =  l 


X  - 


i  =  6 


i=l 


A  further  statistic  computed  in  connection  with  each  X- 
rofile  was  the  correlation  coefficient 


pre 


R  .       =    X 
y/x 


£  "i 

i  =  l 
6 

£  H 

i  =  l 


The  correlation  coefficient  for  each  X-profile  was  general- 
ly in  the  range  .95  to  .99,  which  tended  to  enhance  the 
credibility  of  the  usefulness  of  the  X-profile  technique. 
The  mean  values,  X,  together  with  the  standard  devia- 
tion, O(X),  based  upon  the  25  soundings  on  each  meridian 

are  listed  below: 

Meridian  1     Meridian  2     Meridian  3 


X      3.143 
O(X)      0.799 


3.712 
0.915 


2.931 

1.094 
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C.   PRESSURE-SCALED  ABSORBER  MASSES 

With  the  mixing  ratios  computed  at  each  odd  k-level 
(Figure  1),  it  is  possible  to  compute  the  pressure-scaled 
water  vapor  absorber  masses  in  each  layer.   In  Figure  1, 
these  water  vapor  layer-masses  were  denoted  by  Au,  but 
here  we  develop  notation  which  is  more  representative  of 
the  layers  under  consideration.   For  example,  the  pressure- 
scaled  water  vapor  mass  in  layer  (8,  10)  is  given  by 

q9  Ap    P9     0  72 
•U(8,  10)  =  -_C__^,  •  (2-5) 

where  the  pressure  ratio  raised  to  the  exponent  0.72  (after 
Moller  and  Raschke,  1964)  has  been  employed.   For  the  layer 
(6,  8) ,  an  analog  to  Equation  (2-5)  was  employed  with  q 
replaced  by  q   and  p   by  p  .   A  similar  procedure  was  fol- 
lowed in  the  topmost  layers  (0,  2) , (2,  4)  and  (4,  6) . 
Equation  (2-6)  affords  an  example  for  the  topmost  layer 
(0,  2)  . 


qi  Ap    Pl     0.72 
u(0,  2)  = [I5n7J?]  (2-6) 


In  each  of  the  five  layers  Ap  is  200  mb. 

The  algorithm  for  computing  the  integrated  water  vapor 
mass  above  the  earth's  surface  is  then  as  follows: 

U  (10)  =  0.0 

3)  +  U (8 ,  10) 
)   +  U(6,  8) 

(2-7) 


U(8) 

=    U(10) 

+ 

U (8 ,     10) 

U(6) 

=     U(8) 

+ 

U(6,     8) 

U(4) 

=     U(6) 

+ 

U(4,     6) 

U(2) 

=     U(4) 

+ 

U(2,     4) 

U(l) 

=    U(2) 

+ 

.5 (U(0,      2) ) 

U(0) 

=    U(2) 

+ 

U (0,     2) 
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The  principle  for  the  computation  of  the  carbon  dioxide 
scaled  mass  is  similar  to  that  for  water  vapor.  Here,  how- 
ever, it  is  customary  to  state  the  pressure-scaled  mass  in 

terms  of  the  normal  temperature  and  pressure  (N.T.P.) 

2 
volume  over  a  cm   of  the  earth's  surface.   Thus 


C(8,  10)  =  3.14  x  io"4(_AP  )   ( 


g  p     1013.25 


0.72 


(2-8) 


-4 
In  (2-8)  ,  the  fraction  3.14  x  10    is  the  constant  mixing 


ratio  of  carbon  dioxide  by  volume,  while 


Ap 
9  P. 


is  the  N.T.P, 


thickness  of  the  layer  AP  =  P.  -  P„.   The  quantity  p   is 

108  s 

the  constant  density  in  a  standard  homogeneous  atmosphere 

for  which  T   =  273.16°K  everywhere  and  p   =  1013.25  mb. 
s  ^s 

With  minor  modification,  Equation  (2-8)  becomes 


C(8,  10)  =  3.14  x  10  4 (H) ( 


Ap 


-)  (■ 


1013 .25'  v1013 . 25 


)0'72(2-9) 


RT 


H  = 


7.995  x  10   cm. 


All  other  layer  thicknesses  of  carbon  dioxide  are  similar 
in  form  to  C(8,  10) ,  for  example: 


C  (6,  8)  =  3. 14  x  10~4(H)( 


-)   ( 


1013 .25    1013 .25 


'  )°-72(2-10) 


The  integration  for  the  layer  masses  of  C(8) ,  C(6) ,  C(4) , 
C(2),  C(l)  and  C(0)  follow  exactly  the  procedure  outlined 
for  water  vapor. 
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D.   CLOUD  PARAMETERIZATION 

In  the  calculations  of  the  fractional  cloud  cover 
amounts  in  layers  (4,  6)  and  (8,  9)  the  relative  humidi- 
ties and  thus  the  saturation  vapor  pressure  at  levels 
k  =  5  and  k  =  9  must  be  available.   The  saturation  vapor 
pressures  at  each  level  were  calculated  using  Equations 
(2-11  (a))  and  (2-11 (b)),  and  the  relative  humidity  using 
Equations  (2-12  (a))  and  (2-12 (b)): 

'=<"»'- ii2Hr <2-"<»>> 

qsO00)     =     ■<i2lllS —  <2-ll<b>> 

where  A  =  21.656 
B  =  5418. OK 

R»'5>  =qfrlw  <2-12U" 

R«<5'  -q^oT  (2-12|b" 

The  relative  humidity  values  were  then  used  to  calculate 
the  fractional  cloud  amounts  at  each  layer  using  the  fol- 
lowing formulation,  which  is  similar  to  Smagorinsky ' s  (1960) 

CL(1)  =  2.0(RH(5))  -  0.7  (2-13(a)) 

CL(2)  =  3.33(RH(9))-  2.0  (2-13(b)) 
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III.    SOLAR  RADIATION 


A.   PARTITION  OF  SOLAR  INSOLATION 

In  this  study,  following  Joseph  (1966),  a  solar  con- 
stant of  2.00  ly  min    was  assumed.   This  is  the  value  at 
the  top  of  the  atmosphere;  however,  it  is  assumed  subject 
to  a  four  percent  attenuation  by  oxygen  and  ozone  above 
the  tropopause.   Thus  the  effective  solar  constant  assumed 
at  level  k  =  2  in  this  study  is  1.92  ly  min 

The  effective  solar  insolation  at  the  tropopause  was 
computed  from 

F(2)  =  S[^-]~2  cos  Z  (3-1) 

R 
m 


where 

5  =  effective  solar  constant 

Cos  Z  =  cosine  of  the  zenith  angle 

R/R    =  ratio  of  the  actual  earth-sun  distance 
to  the  mean  earth-sun  distance. 

The  ratio  R/R   is  listed  as  a  function  of  the  Julian  date 
m 

in  the  Smithsonian  Meteorological  Tables  (List,  1958). 
In  this  study  based  upon  data  for  25  April  at  0000GMT  a 

value  of  R/R   =  1.00606  was  used.   The  cosine  of  the 

m 

zenith  angle  was  calculated  by 

Cos  Z  =  Sin4»  Sin<5  +  Cos<J>  Cos6  Cos  h         (3-2) 

where  $    -    Latitude 

6  =  Solar  declination  for  April  25th  =  12°56' 
h  =  hour  angle 

22 


The  value  of  the  declination  was  taken  from  Table  169 
(List,  1958),  and  the  value  of  Sin((>  was  calculated  using 
two  different  formulas  depending  on  the  data-line  used 
for  the  computations.  These  values  of  Sin<J>  along  lines 
one,  two,  and  three  are  calculated  in  terms  of  the  FNWC 
map  coordinates  (I, J)  as: 


Lines  1,3     Sin<f>  = 


Line  2 


Sin<f>  = 


973.752  -  2(32-1)' 

973.752  +  2(32-1)' 

9  7  3.752  -  (32-1) 2 

973.752  +  (32-I)2 


(3-3a) 


(3-3b) 


I  =  1,  ,  25 

Here  I  is  the  abscissa  distance  on  the  FNWC  grid,  Figure  2. 
The  soundings  were  all  taken  at  OOOOGMT;  therefore,  the 
value  for  the  hour  angle  h,  can  easily  be  calculated.   For 
example,  in  this  experiment  solar  noon  existed  at  the  180th 
meridian  and  the  hour  angles  h  =  55°,  10°  and  35°  occurred 
along  lines  one,  two,  and  three,  respectively. 

Following  Joseph  (1966)  a  computationally  simple  parti- 
tion of  the  solar  insolation  was  used  for  this  study. 
Joseph  divided  the  insolation  at  level  2,  F(2) ,  into  two 
parts,  one  of  which  was  subject  to  water  vapor  absorption 
but  not  to  Rayleigh  scattering.   This  portion  is  given  by 


F (A)  =  0. 3  4  9F  (2)  ; 


(3-4) 
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the  remaining  portion  of  the  insolation, 


F(S)  =  0.651F (2) , 


(3-5) 


is  subject  only  to  Rayleigh  scattering  (but  not  to  water 
vapor  absorption)  in  a  clear  moist  atmosphere.   The 
presence  of  two  cloud  decks  in  this  study  introduces 
cloud-reflectivity  into  both  the  F(A)  and  F(S)  solar 
energy  portions. 

The  F(S)  energy  is  composed  of  all  wavelengths  A  < 
0.9  Urn  since  here  water-vapor  absorption  is  negligible. 
The  F (A)  portion  consists  of  those  wavelengths  A  >  0.9  urn 
where  absorption  by  water-vapor  and  carbon  dioxide  bands 
in  the  near  infra-red  is  the  domine.nt  attentuation  process 
in  clear  air. 

B.    DISPOSITION  OF  F(S)  INSOLATION 

Joseph  (1966)  found  that  empirical  values  of  Rayleigh 
clear-sky  albedo,  (after  Coulson,  1959),  were  best-fit 
using  the  least  squares  technique  by  the  following  formula 


a(R)  =  .085  +  .25074  [log,„  (—  Sec  Z)] 

10   pQ 

p   =  1013.25 


(3-6) 


Here  T/p   =  1  in  view  of  the  fact  that  mean  sea  level 
occurs  near  1000  mb .   The  value  Sec  Z  is  equal  to  (Cos  Z) 
where  Cos  Z  is  calculated  from  (3-2)  . 

Another  reflective-type  parameter  which  must  be  con- 
sidered in  this  study  is  the  surface  albedo  a(G) .   Over 


-1 
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the  land  areas  of  the  Northern  Hemisphere  a  constant 
average  surface  albedo  of  0.14  is  utilized.   However,  the 
value  over  the  oceans  depends  on  cos  Z.   Since  in  this 
study  all  of  the  soundings  were  over  the  ocean,  the 
equation 

<X(G)  =  max  {.06,  .06  +  .54  (.7  -  cos  Z)  }   (3-7) 

was  used  to  calculate  a(G),  (after  Gates  et  alT  (1971)). 

Since  a  two-layer  cloud  model  was  used,  four  separate 
cloud  cases  are  discussed  below. 

1 .   Clear  Sky  Case 

In  this  case  the  F(S)  insolation  is  subject  only 
to  Rayleigh  sky  reflectivity,  a(R),  and  the  earth-surface 
reflectivity,  a(G),  of  Equations  (3-6)  and  (3-7) „   Con- 
sidering the  possibility  of  multiple  reflections  between 
the  earth  and  atmosphere,  each  of  which  diverts  downward 
the  fraction  a(R)  of  the  earth's  surface  reflectance  a(G), 
the  clear  sky  insolation  at  the  earth's  surface  after 
scattering  is  given  by 

IS10(0,0)  =  F  (S)  [l-a(R) ]  [l-a(G) ]/ [l-a (R)a(G) ]  (3-8) 

In  the  next  three  cases,  where  clouds  were  present, 
the  cloud  reflectivity  values  suggested  by  C.  D.  Rodgers 
(1967)  were  used.   Thus  for  the  higher  level  clouds  (between 
k  =  4  and  6) ,  R(l)=  0.54  and  for  the  lower  level  clouds  (be- 
tween k  =  8  and  9),  R(2)  is  0.66. 
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2.   Cloudy  Sky  Cases 

We  consider  first  the  fraction  of  the  grid  area 
having  clouds  at  both  levels.   The  F(S)  insolation  ab- 
sorbed by  the  ground  at  such  a  grid  point  is  given  by 


IS10(1,1)  =  F(S)  [  (l-R(l)  )  (1-R(2)  )  (l-Ct(G)  )] 
*[l-((R(l)R(2)+R(2)Cl(G)+R(l)a(G)) 


(3-9) 


+  2(R(l)R(2)a(G)  )  ] 


-1 


(after  Arakawa , 1972 ) .   The  special  symbols  (1,1)  denote 
CL(1)  =  CL(2)  =  1.0,  which  will  temporarily  be  assumed. 

The  other  cloud-cover  variations  to  be  considered 
are  those  of  CL(1)  =  1.0,  CL(2)  =  0.0  and  CL ( 1 )  =  0.0, 
CL(2)  =  1.0,  respectively.   In  the  former  case,  R(2)  = 
0.0  in  (3-9)  so  that 

IS10(1,0)  =  F(S)  (l-R(l)  )  (l-a(G)  )/(l-R(l)a(G)  )    (3-10) 

In  the  latter  case  with  clouds  at  the  lower  layer  only  it 
follows  that 

IS10(0,1)  -  F(S)  (1-R(2)  )  (l-a(G)  )/(l-R(2)a(G)  )    (3-11) 

Equations  (3-8),  (3-9),   (3-10),  and  (3-11)  were  used  to 
calculate  the  absorbed  F(S)  insolation  at  the  earth's  sur- 
face after  only  reflective- type  processes  are  considered 
to  take  place.   However,  in  each  case  the  sky  was  con- 
sidered to  be  either  overcast  or  clear  and  an  appropriate 
grid-area  weighting  scheme  was  necessary  to  correct  to  the 
actual  cloud  conditions  CL ( 1 ) ,  CL ( 2 ) . 
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3 .   The  Weighted  F(S)  Insolation 

The  cloud  combinations  denoted  by  (0,0),  (1,1), 
(1,0),  and  (0,1),  respectively,  actually  occur  with  the 
areal  coverage  (or  weights)  at  a  grid  point  given  by  the 
following  functions: 


W(0,0)  =  (l-CL(l) ) (1-CL(2) ) 

W(l,l)  =  (CL(1)  )  (CL(2)  ) 

W(1,0)  =  (CL(1)  )  (1-CL(2)  ) 

W(0,1)  =  (l-CL(l)) (CL(2) ) 


(3-12) 


at  each  grid  point. 

When  these  weights  were  applied  multiplicatively 
in  sequence  with  the  right  sides  of  Equations  (3-8),  (3-9), 
(3-10),  and  (3-11),  then  the  composite  surface-absorbed 
insolation,  IS10  resulted. 


IS10  =  IS10  (0,0) W(0,0)+  IS10 (1,1) W(l#l) 
+  IS10 (1,0) W(1,0)+  IS10 (0, 1) W(0,1) 


(3-13) 


Note  that  by  subtraction  of  IS10  from  F(S)  we  also 
arrive  at  that  portion  of  the  F(S)  insolation  reflected  to 
space . 


C.    DISPOSITION  OF  F(A)  INSOLATION 

In  the  following  sections,  the  fraction  of  the  insola- 
tion which  is  absorbed  by  atmospheric  water-vapor  and 
carbon  dioxide  is  discussed. 

1 .   The  Clear  Sky  Case 

In  this  connection  the  Manabe-Mol ler  absorptivity 
function 


a(2,2k)=  . 271 [U (2 , 2k) Sec  Z] 


0.  303 


(3-14) 
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Figure  3.   Schematic  representation  of  F(A)   insolation  dis- 
position in  the  clear  sky  case. 
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was  used.   The  form  of  Equation  (3-14)  indicates  that  the 
absorptivity  a_  is  applied  over  the  pressure-scaled  water- 
vapor  mass  between  levels  2  and  2k  (see  Figure  1)  along 
the  zenith  slant-path  angle  Z.   Therefore,  assuming  negli- 
gible water-vapor  mass  above  level  2,  the  absorbed  solar 
insolation  in  the  layer  (2,6)  is 


A(2,6)   =  0.271F(A) [U(2,6) Sec  Z] 


0.303 


(3-15) 


Similarly,  the  absorbed  insolation  in  the  layer  (2,10)  is 


A(2,10)  =  0.271F (A) [U(2,10)Sec  Z] 


0.303 


(3-16) 


and  the  absorbed  insolation  in  the  layer  (6,10)  is  given 
(see  Figure  3)  by 


A(6,10)  =  A (2 ,10)  -  A (2,6) 


(3-17) 


The  direct  transmission  of  F(A)  insolation  to  the 
earth's  surface  was  found  by  subtracting  A(2,10)  from  F(A). 
This  direct  transmission  is  partially  reflected  by  the 
earth's  surface  (the  value  computed  for  a(G)  in  Equation 
(3-7)  is  applicable  here)  and  the  remainder,  absorbed  at 
the  earth's  surface,  is  given  by 

IA10(0,0)=F(A){l-.271[U(2,10)Sec  Z]°*303}(l-a(G)) 

(3-18) 
In  the  following  subsections  dealing  with  partially 
cloudy  layers,  it  is  necessary  to  adopt  representative 
cloud  reflectivities.   For  this  purpose,  values  suggested 
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by  C.  D.  Rodgers  (1967)  were  used.   These  reflectivity 
values  for  the  F(A)  wavelengths  differ  from  those  listed 
by  Rodgers  for  the  F(S)  wavelengths.   We  followed  the 
procedure  of  considering  overcast  cloud  conditions  when- 
ever clouds  occur  and  subsequently  applying  the  weight- 
ing factors  as  was  done  in  Section  III(B). 
2 .   Overcast  in  Both  Cloud  Layers 

Here,  with  clouds  at  both  levels  the  cloud  reflec- 
ivity  values  RA(1)  =  .46  and  RA ( 2 )  =  .50  were  used.   In 
addition  Rodgers  suggests  cloud  absorptivities  of  A(l)  = 
.20  and  A(2)  =  .30  (see  Figure  4).   With  the  following 
set  of  formulas,  the  downward  directed  insolations  FA44-, 
etc.,  were  calculated  at  each  required  level  in  order  to 
arrive  at  a  value  for  the  insolation  abosrbed  by  the 
atmospheric  layers,  and  finally  by  the  earth's  surface. 


0.303 


A(2,  4)'  =  F  (A)  . 271 [U (2,4) Sec  Z] 

FA44-       =  F  (A)  [1-.  271  (U(2,4)  Sec  Z)°*3°3] 


A(2,  4)"  =  F (A)RA(l)  . 271  [U  (2,4) Sec  Z] 
A ( 2 ,  4 )   =A(2,  4 ) '  +  A ( 2 ,  4)" 
A(4,  6)    =  FA4iA(l) 

=  FA44-  (1-RA  (1)-A  (1)  ) 


0.303 


FA64- 
FA8I 


(3-19a) 

(3-19b) 

(3-19c) 

(3-19d) 

(3-19e) 

(3-19f ) 

FA64-(TD(6,  8)  )  [1-RA  ( 1  )  RA  (  2  )  (TD  (6,  8))2]"1 

(3-19g) 


A(6,  8)  =  FA64-  -  FA8I 

A  (8,  9  )  =  FA8IA  (2) 

FA94-  =  FA84-  (1-RA  (2) -A  (2)  ) 

FA104- 


(3-19h) 

(3-19i) 

(3-19J) 
2.  -1 


=  FA9+TD(9,  10)  [1-CX  (G)  RA  (2)  (TD  (9,  10))   ] 

(3-19k) 
A(9,  10)  =  FA94-  -  FA10+  (3-19P,) 
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Figure  4.   Schematic  representation  of  F(A)  insolation  dis- 
position in  the  case  of  two  overcast  layers. 
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where 


TD(6,  8)  =  1  -  .271 [U(6,8) 5/3] 


0.303 


(3-20) 


TD(9,  10)=  1  -  .271  [D(9,10) 5/3] 


0.303 


In  the  preceding  equations  A(2,  4),  A(4,  6), 
A(6,  8),  A(8,9),  and  A(9,10)  represent  the  insolation  ab- 
sorbed in  the  indicated  layers.   A  superscript  prime  fol- 
lowing A(2,4)  indicates  absorption  (by  use  of  Equation  3-14) 
of  the  downward  directed  beam,  while  the  double  prime  indi- 
cates absorption  of  the  beam  specularly  reflected  upward 
from  the  cloud  top.   The  functions  of  Equation  (3-20)  repre- 
sent transmissivity  after  diffuse  passage  through  a  cloud 
or  sub-cloud  layer,  which  causes  the  effective  replacement 


of  Sec  Z  in  Equation  (3-14)  by  Sec  Z  =  5/3  (Katayama,  1966). 
The  symbols  FA64- ,  FA84-,  FA9+,  and  FA10  +  indicate  the  down- 
ward streams  of  insolation  passing  through  the  indicated 
levels.   Note  that  the  effect  of  multiple  reflections  be- 
tween clouds  is  incorporated  in  the  denominator  of  FA8+, 
while  that  between  the  lower  cloud  and  the  ground  has  a 
corresponding  term  in  the  denominator  of  FA10+,  involving 
the  reflectivities  a(G)  and  RA ( 2 ) .   The  effect  of  the 
multiple  reflectances  as  modeled  here  is  to  enhance  the 
insolation  at  the  lowermost  reflecting  boundary  of  the  two 
under  consideration. 

Finally,  note  that  the  flux  absorbed  by  the  earth's 
surface  under  the  conditions  of  an  overcast  layer  at  both 
levels  is  given  by 


IA10 (1,1)  =  FA10 (l-a(G) ) 


(3-21) 
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3 .   Disposition  of  F(A)  Insolation  with  an  Upper 
Overcast 

In  this  case  with  only  an  upper  overcast  layer, 

the  values  of  RA ( 1 )  and  A(l)  used  are  the  same  as  those 

of  the  preceding  section.   However,  a  somewhat  simpler 

set  of  disposition  equations  result: 

A(2,  4)=  F(A) {.271[U(2,  4)Sec  Z]°*3°3 

+[1  -0271(U(2,  4)Sec  Z) ° * 3°3 ] RA ( 1 ) }  (3-22a) 
A(4,  6)=  F(A)A(1) [1-.271(U(2,  4)Sec  Z)°'3°3]  (3-22b) 
FA6i  =  F  (A)  [1-.271  (U(2,  4)Sec  Z)  *  3  °3  ]'  [  1-RA  ( 1 )  -A  ( 1 )  ] 


(3-22c) 


FA10  4-  =  FA6+  [1-.271 (U(6,  10)5/3)°*3°3] 


*  [l-a(G)RA(l)TD (6,  10)2]"1  (3-22d) 

A(6,  10)=  FA64-  -  FA10  +  (3-22e) 


Here  the  variables  are  defined  similarly  to  those  of  the 
preceding  section.    (Note  that  A(2,  4) '  and  A  (2,  4)"  have 
been  combined  in  Equation  (3-22)  to  give  A(2,  4)  and  the 
definition  of  TD(6,  10)  is  similar  in  principle  to  that 
of  TD(6,  8)).   The  absorbed  insolation  at  the  ground  is 
given  by 

IA10(1,  0)  =  FA104-   (1  -  a(G))  (3-23) 

(see  Figure  5). 

4 .   Disposition  of  F(A)  Insolation  with  a  Low  Overcast 
The  set  of  equations  which  describe  the  behavior 
of  the  insolation  with  a  low  overcast  deck  are  given 
below  (see  Figure  6) . 
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Figure  5.   Schematic  representation  of  F(A)  insolation  dis- 
position with  an  upper  overcast  layer  only. 
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Figure  6.   Schematic  representation  of  F(A)  insolation  dis- 
position with  a  lower  overcast  layer  only. 
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0.303 


0.303 


A(2,  6)'  =  .271F (A) [U(2,  6)Sec  Z] 

A(6,  8)'  =  .271F(A) [U(2,  8)Sec  Z) 

-(U(2,  6)Sec  Z)°*303] 

FA8+  =  F (A) [1-.271 (U(2,  8)Sec  Z)°*3°3] 


FA9  +  =  FA8+ [1-RA(2)-A(2) ] 

FA8t  =  FA8iRA(2) 

A(6,  8)"  =  FA8+ [ .271 (U (6,  8)Sec  Z) 

A  ( 6  ,  8  )   =A(6,  8  )  '  +  A  ( 6  ,  8)" 


0.303 


A(2,  6)"  =  .271FA8T [ (U(2,  8)Sec  Z) 


0.303 


-(U(6,  8)Sec  Z)0,303] 


A(2,  6)  =  A(2,  6)'  +  A(2,  6)" 
A(8,  9)  =  FA84-  A(2) 


(3-24a) 

(3-24b) 
(3-24c) 
(3-24d) 
(3-24e) 
(3-24f ) 
(3-24g) 

(3-24h) 
(3-24i) 
(3-24J) 


FA104    =  FA9+TD(9,  10 ) [ 1-a (G) RA ( 2 ) (TD ( 9 ,  10))2]  1 


A(9,  10)=  FA9+  -  FA10* 


(3-24k) 
(3-24£) 


These  variables  also  have  similar  definitions  to  those  in 
Section  III-C-2.  As  before  the  absorption  by  the  earth's 
surface  is  given  by 


IA10(0,  1)  =  FA104-  (1  -  a(G)) 


(3-25) 


5 .   Weighted  F(A)  Surface  Insolation 

As  previously  described  in  the  section  on  the 
F(S)  insolation,  the  values  of  IA10  calculated  for  the 
cases  (0,  0),   (1,  1),   (1,  0)  and  (0,  1)  must  be  appro- 
priately weighted  by  the  values  described  in  Equation  (3-12) 
in  order  to  arrive  at  the  composite  F(A)  surface-absorbed 
insolation . 
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IA10  =  IA10(0,  0)W(0,  0)+  IA10(1,  1)W(1,  1) 

+  IA10(1,  0)W(1,  0)+  IA10(0,  1)W(0,  1)      (3-26) 

6 .   Atmospheric  Layer  Absorption  of  F (A) 

In  each  of  the  four  cases  discussed,  the  F(A)  in- 
solation was  absorbed  as  it  passed  through  the  key  atmos- 
pheric layers  (2,  6)  and  (6,  10).   In  the  clear-sky  case 
these  absorptions  are  given  by  Equations  (3-15)  and 
(3-17).   However,  in  the  calculations  of  the  cloudy-sky 
cases  the  various  contributions  have  merely  been  added 
together  to  give  the  total  insolation  absorbed  in  the 
layer  (2,  6)  and  then  in  the  layer  (6,  10) ,  even  though 
cloudiness  affects  the  radiation  physics  differently  within 
the  cloud  and  subcloud  layers  as  Figures  3,  4,  5,  and  6 
indicate . 

a.   Overcast  in  Both  Cloud  Layers 

In  this  case  the  insolation  absorbed  in  the 
two  layers  was  calculated  with  values  from  the  appropriate 
equations  from  (3-19a,  ...,  3-19£) 


A(2,  6)      =  A(2,  4)'  +  A  (  2  ,  4)"  +  A(4,  6) 

-L  »  -L 

A(6,  10)     =  A(6,  8)    +  A(8,  9)    +  A(9,  10) 
1»1 


(3-27) 


where  the  subscript  (1,  1)  on  A(6,  10)  indicates  the  cloud 
layer  combinations  which  will  subsequently  be  weighted, 
b.   Upper  Overcast  Only 

For  the  high  level  overcast  only  the  insola- 
tion absorbed  was  calculated  using  the  various  equations 
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(3-22a,  . . . ,  3-22e) 


A(2,  6^  Q   =  A(2,  4)  +  A(4,  6) 

A(6,  10)      =  FA64-  -  FA10  + 
1,0 


(3-28) 


c.   Lower  Overcast  Only 

In  the  low  level  overcast  cas.e  the  value  of 
the  insolation  absorbed  in  the  layers  ( 2  „  6)  and  (6,  10) 
was  calculated  using  equations  of  (3-24a„...,  3-245,) 


A  (  2  ,  6  )       =  A(2,6)'+A(2,  6)" 

u  ,  i 


A(6,  10) 


0,1 


(3-29) 
A(6,  8)'  +  A(6,  8)"  +  A(  8  ,  9)  +  A(9,  10) 


d.   Composite  Weighted  Absorption  by  Layers 

In  this  case,  as  in  the  previous  subsection 
on  weighted  F(A)  surface  insolation,  the  values  computed 
in  each  of  the  four  cases  must  be  appropriately  weighted 
by  the  functions  described  in  Equation  (3—12)  in  order  to 
arrive  at  a  composite  value  of  the  F (A)  insolation  absorbed 
in  the  atmospheric  layers  (2,  6)  and  (6,  10). 


A(2,  6)=A(2,  6)Q  Q  W(0,  0)+A(2,  6)      W(l,  1) 

+A(2,  6)1    Q  W(l,  0)+A(2,  6)Q  1    W(0,  1)         (3-30) 

A(6,  10)=A(6,  10)    QW(0f  0)+A(6,  10)    jW«l,  1) 

+  A(6,  10)    0wd,  0)+A(6,  10    W(0,  1)  (3-31) 

e „   The  Absorptivity  (ABA)  by  Layers 

In  this  study  it  was  convenient  to  deal  with 
the  fractional  absorptivity,  or  simply  absorptivity,  as 


38 


well  as  the  actual  insolation  energies  absorbed  in  each 
layer,  given  by  Equations  (3-32  and  3-33)  and  by  Equation 
(3-26)  at  the  ground.   In  the  calculation  of  the  absorp- 
tivity, the  total  incoming  insolation  at  the  top  of  the 
atmosphere  must  be  known.   This  value  was  calculated  us- 
ing Equation  (3-32) 

FADJ  =  2.00(R/R  )~   Cos  Z  (3-32) 

m 

Table  II  lists  the  instantaneous  values  of  FADJ  at  .each 
gridpoint  along  the  three  data  lines.   The  calculation  of 
the  f ull-tropospheric  absorptivity  then  follows  from 

,„„    A(2,  6)  +  A(6,-  10) 

ABA  =  ¥ZET-* (3"33) 

which  is  the  ratio  of  the  insolation  absorbed  in  the  tropo- 
sphere to  the  extraterrestrial  insolation  at  the  same  time 
for  each  gridpoint. 

D.    REFLECTION  OF  F(A)  AND  F(S)  ENERGY 

In  the  calculations  of  the  albedo  and  the  fractional 
absorptivity  of  the  ground  both  the  F(A)  and  F(S)  frac- 
tions of  the  solar  insolation  had  to  be  considered. 

1 .   Albedo 

That  portion  of  the  F(A)  insolation  which  was  re- 
flected back  to  space  is  given  by 

REFA  =  [F(A)  -  A(2,  6)  -  A(6,  10)  -  IA10]         (3-34) 
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TABLE  II.   Total  extraterrestrial  insolation  (FADJ)  at 
each  grid-point  of  the  three  data  lines  (ly  min~l). 


Line  1 


Line  2 


Line  3 


.899 
.923 
.947 
.970 
.995 
1.018 
1.042 


1 

1, 

1, 

1. 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1 


064 
085 
106 
125 
142 
158 
171 
180 
187 
190 
187 
181 
168 
149 
122 
088 
048 
998 


1.899 
1.912 
923 
933 
939 
946 
947 
947 
942 
934 
923 
907 
888 
862 
834 
798 
757 
710 
658 
598 
534 
462 
384 
301 
210 


1.346 
1.374 
1.403 
1.430 
1.458 
485 
510 
534 
558 
579 
597 
613 
626 
634 
638 
637 
1.629 
1.616 
593 
565 
526 
478 
419 
350 
270 
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where  the  calculation  equations  of  A(2,  6> )  ,  A(6,  10),  and 
IA10  have  been  described  in  Equations  (3—30),  (3-31),  and 
(3-26),  respectively. 

The  F(S)  insolation  reflected  was  simpler  to  calcu- 
late because  the  F(S)  wavelengths  do  not  undergo  atmospheric 
absorption  in  the  passage  through  the  layers.   The  value 
of  the  reflected  F(S)  insolation  is  given  by 


REFS  =  F(S)  -  IS10 


(3-35) 


where  IS10  was  calculated  by  Equation  (3—13). 

These  values  of  REFA  and  REFS  were  then  combined 
to  give  a  value  for  REF,  which  is  the  total  reflected  energy 
passing  outward  through  the  tropopause  (k.=  2) 


REF  =  REFS  +  REFA 


(3-36) 


The  albedo  is  defined  as  the  ratio  of  the  total  re- 
flected energy  to  the  total  incoming  energy.   This  ratio 
was  calculated  by 


ALB  = 


REF 
FADJ 


(3-37) 


where  FADJ  is  defined  in  Equation  (3-32)  .   The  values  of 
REFA  and  REFS  were  computed  at  each  grid-point.   The  re- 
sulting REF  values  are  shown  in  cross-section  form  in 
Figures  11(a)  and  11(b).   Also,  values  of  the  albedo  were 
computed.   The  global  mean  albedo  over  the  75  soundings 
turned  out  to  be  .412,  based  upon  the  cloud-cover  para- 
meterization for  CL(1),  CL(2). 
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2.  Composite  Absorptivity  (ABG)  by  Earth 

The  total  insolation  absorbed  at  the  earth  was 
the  sum  of  the  weighted  values  of  the  F(S)  and  the  F (A) 
portions  of  the  incoming  solar  insolation. 

Q(E)  =  IA10  +  IS10  (3-38) 

Here  Q(E)  is  the  insolation  absorbed  by  the  earth's  sur- 
face at  each  grid-point  sounding. 

Using  Q(E),  a  fractional  absorptivity  for  the 
earth's  surface  was  calculated.   The  equation  used  for 
the  calculation  was 

where  the  calculations  of  Q(E)  and  FADJ  iiave  been  previously 
discussed „ 

3 .  Check  on  Computations 

The  actual  computational  process  used  in  this  model 
was  checked  against  the  condition  that  the  summation  of  the 
fractional  values  ALB,  ABA,  and  ABG  at  each  grid-point 
equals  .96.   This  must  be  expected  since  we  allowed  for 
attenuation  of  the  solar  insolation  in  the  amount  of  four 
percent  as  it  passed  through  the  stratosphere.   This  check 
computation  was  satisfied  at  each  grid-point. 

E.    STATISTICAL  COMPUTATIONS 
1 .   Clear  Sky  Cases 

In  order  to  verify  the  preceding  calculations, 
statistical  computations  of  several  key  values  were  obtained 
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using  linear  regression  techniques  from  the  BMD  set  of 
statistical  programs  (Dixon,  1973).   The  predictand  of 
the  first  computation  was  ALB  (0,0, Z),  which  is  the 
albedo  in  the  clear  sky  case.   The  predictor  was  the  se- 
cant of  the  zenith  angle  (Sec  Z).   The  best-fit  equation 
which  resulted  was 


ALB(0,0,Z)  =  (-0.02726  +  0.11387  Sec  Z) 


(3-40) 


with  the  regression  coefficient  R  =  .979-   Equation  (3-40) 
summarizes  how  the  radiative  model  linearizes  the  effect  of 
Sec  Z  in  the  clear  sky  case.   The  high  R— value  confirms 
the  known  result  that  Sec  Z  is  a  good  predictor  of  the 
albedo  (with  clear  skies),  which  is  to  be  expected  over 
an  oceanic  area.   The  mean  albedo  in  the  clear  sky  case, 
which  corresponds  to  the  average  secant  of  the  zenith 
angle,  is  as  follows 


ALB  (0,  0)  =  .135 


Sec  Z  =  1.428 


Two  other  regressions  that  were  performed  using 
the  linear  regression  technique  were  ABA  (0,0, M)  and 
ABG  (0,0, M)  against  the  water-vapor  mass  path  length  M, 
where  M  is  given  by 


M  =  (U  Sec  Z)  ' [log     (U  Sec  Z)  /2] 


This  parameterization  of  the  water  vapor  mass  is  essential- 
ly that  of  Hanson  (1971),  who  used  a  similar  absorber-mass 
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predictor  in  his  empirical  specification  of  ABA  for  both 
clear  and  partly  cloudy  skies.   The  resulting  equations 
obtained  with  clear  skies  in  this  study  for  ABA  (0,0, M) 
and  ABG  (0,0, M)  were 

ABA  (0,0, M)  =  .096  +  .031(U  Sec  Z)1/2[logj   (U  Sec  Z)1'2] 
R  =  .984  (3-41) 

1/2  1/2 

ABG  (0,0, M)  =  .800  -  .035(U  Sec  Z)  '      [logj   (U  Sec  Z)  '   ] 

R  =  .956  (3-42) 


The  means  and  standard  deviations  of  ABA  (JO ,  0  ,M)  ,  ABG(0,0,M), 
and  M,  based  on  all  of  the  original  75  soundings ,  are 
given  below 


ABA(0,0,M)  =  .137  O    =     .029 


ABG(0,0,M)  =  .687  O    =     .034 

M   =  1.342 (gm  cm"2)1^2     a  =  .93  4  . 


2 .   Relationship  between  Albedo  in  the  Cloudy  and 
Clear  Sky  Cases 

As  discussed  previously,  it  was  pxossible  to  compute 
the  value  of  ALB  in  the  clear  sky  case  amcd  also  in  the  com- 
posite cloudy  sky  case.   In  order  to  stuctry  the  relationship 
between  these  two  cases  a  ratio  ALB  [CL  ( 1 )' ..,  CL  (  2  )  ] /ALB  (0  ,  0) 
for  the  same  zenith  angle  at  grid  point  (.  K  ,  J )  was  formed. 
This  value  was  used  as  the  predictand  and!  the  total  opaque 
c loud  cover  (CL) ,  specified  by  the  model  was  used  as  the 
predictor.   Here  CL  is  calculated  using  Equation  (3-43), 
which  approximates  the  effective  cloud-cover  by  the  two 
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layers  of  clouds: 

CL  -  CL(1)  +  CL(2)  -  CL(1)(CL(2)     (3-43) 

The  regression  equation  which  resulted  was 

ALB[CL(1)  ,CL(2)  ]  =  ALB(0,0)(1  +  3.931CL  -  . 87 82CL2 )  (  3-44 ) 
R  =  0.956 

The  means  and  standard  deviations  of  the  values  in  Equa- 
tion (3-44)  are 


ALB  [CL (1)  ,CL (2) ]  =  .412  a    =  .431 


ALB(0,0)  =  .135  a  =  .141 

CL  =  .680  a  =  „742 

Note  that  the  mean  cloudy  sky  albedo  generated  by 
the  model  is  3.211  times  that  of  the  mean  clear  sky  albedo. 
The  high  value  of  the  correlation  coefficient  (R  =  .956) 
indicates  reliability  of  the  indicated  dependence  upon  CL. 
This  global  mean  albedo  of  .412  is  somewhat  higher  than 
that  recently  reported  [ALB  =  .300]  by  Raschke  et  al  (1973) 
However,  the  mean  value  reported  here  does  not  have  an 
input  from  the  polar  latitudes.   Raschke,  et  al  suggest 
that  perhaps  the  reduced  cloud  amounts  (CL)  recently  re- 
ported in  low  latitudes  would  account  for  the  reduction 
compared  to  that  computed  here  by  the  parameterization 
Equation  (2-13),  after  Smagorinsky. 


45 


3 .   Statistical  Relationship  between  Ground  Insolation 
in  the  Cloudy  and  Clear  Sky  Cases 

In  this  section  as  in  the  preceding  one  a  ratio 

Y  =  ABG (CL) /ABG (0 ,  0)  was  formed  for  each  sounding  case 

to  determine  the  best-fit  regression  formula  using  CL  as 

the  independent  variable.   The  formula  which  resulted  was 

Y  5  ABG  (CD/ABG  (0, 0)  =  (1  -  .66CL),   R  =  .999   (3-45) 

The  ratio  of  ABG (CL) /ABG (0 ,  0)  is  essentially  the 
corresponding  ratio  of  the  transmissions  received  by  a 
pyrheliometer .   Quinn  (1970)  has  summarized  the  empirical 
solar  tr ansmissivity  formulas  by  many  investigators  and 
shows  that  some  of  the  most  reliable  estimates  are  very 
nearly  as  formulated  by  Equation  (3-45).  For  example, 
Budyko  (1956)  for  the  tropics  lists 

ABG (CL) /ABG (0,  0)  =  1  -  .655CL 
and  Kimball  (1928)  lists  the  corresponding  ratio  as 

ABG (CL) /ABG (0,  0)  =  1  -  .71C 
The  mean  statistics  over  the  sample  soundings  gave 


ABG(CL)  =  .379 


ABG(0,  0)  =  .688 
cT   =  .680 
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4 .   Absorptivities  in  Air 

Here  the  computed  values  of  the  ratio  ABA(CL)/ 
ABA(0,  0)  were  investigated  by  least  squares  with  the 
result 

ABA (CL) /ABA (0,  0)  =  1  +  .37863CL,  R  =  .781    (3-46) 

with  mean  statistics  for  the  75-case  population 


ABA(CL)     =  .169 


ABA(0,  0)  =  .137 
CL   =  .680 

This  indicates  that  the  atmospheric  absorptivity  is  some- 
what enhanced  with  the  cloud-model  solar  dispositions  in- 
troduced.  An  increasing  solar  absorptivity  with  increasing 
cloud  cover  is  in  general  agreement  with  Plante's  results 
(1972)  in  a  one-layer  cloud  model. 
5 .   General  Conclusions 

The  statistical  relationships  developed  from  the 
radiation  model  show  consistency  with  results  reported  by 
other  researchers.   For  example,  Equation  (3-41)  gives 
ABA  which  closely  fits  Hanson's  (1971)  empirical  results, 
and  Equation  (3-45)  gives  a  ratio  of  transmittance  which 
closely  agrees  with  the  empirical  results  of  Budyko  (1956) 
and  Kimball  (1928).   Finally,  the  mean  value  of  ABA  in 


clear  skies  (ABA  =.137)  corresponds  closely  to  the  tropo- 
spheric  mean  absorptivity  of  15  percent  reported  by  London 
(1957)  . 
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IV.   TERRESTRIAL  RADIATION 

A.   THEORETICAL  AND  EMPIRICAL  BASIS 

Sasamori  (1968)  developed  empirical  formulas  for 
computing  emissivities  for  long-wave  flux  calculations 
associated  with  the  NCAR  General  Circulation  Model. 
These  empirical  emissivity  formulas  are  matched  to  the 
values  developed  by  Yamamoto  (1952)  from  radiative  transfer 
theory,  which  have  also  been  built  into  his  well-known 
Radiation  Chart.   This  chart  has  proved  to  be  quite  accu- 
rate for  computational  purposes,  and  is  used  here  as  a 
schematic  guide  for  integration  of  the  radiative  transfers 
through  the  various  layers  of  interest  in  the  soundings 
discussed  in  Section  II0 

For  utilization  in  the  FNWC  heating  package  the  neces- 
sary emissivity  formulas  of  Sasamori  are  those  which  will 
lead  to  expressions  for 

F     =  net  IR  flux  at  earth 

* 
F.    =  net  IR  flux  at  level    k  =  6 

D 

* 

F„    =  net  IR  flux  at  level    k  =  2 
2 

*     *         *     * 
The  differences  F   -  F    and  F   -  F   will  then  give  the  IR 

flux-divergence  (i.e.,  cooling  rates)  in  the  layers  (6,  10) 

and  (2,  6),  respectively.   A  workable  scheme  for  making 

such  computations  is  presented  below  for  various  combinations 

of  cloud  cover  [CL(1),  CL ( 2 ) ] . 
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For  example,  in  order  to  compute 

B, 


10  =  Bio  "/ 


10 


£   (U,  C,  T)  dB 
wc 


(4-1) 


B  =  0 


one  must  have  a  good  representation  of  the  emissivity  £ 


wc 


as  a  function  of  both  water  vapor  and  CO_  absorber  masses 
along  the  sounding,  that  is  along  the  path  of  integration. 
Figure  7(a)  illustrates  the  concept  involved  at  level 
k  =  10,  with  clear  skies,  where  the  hatched  area  represents 
the  downward  flux  at  k  =  10.   Here  the  radiative  sounding 
is  converted  into  emissivities  as  indicated  in  the  figure. 
The  hatched  area  denoted  by 


B 


10 


/ 


B=0 


£    dB 

wc 


is  then  the  key  quantity  to  be  evaluated.   For  this  integral 
evaluation,  Sasamori  (1968)  proposes  emissivity  formula- 
tions of  the  type 

£    (8,  10)  =   i  .240  login[0(8,  10)]  +  .622V 

WC  [  1 U  ) 


+     .07262  < (1 


*  [log1Q  C(8,  10)  +  1.064]  \ 


62556[U(8,  10)  +  .0286]*26) 

(4-2) 


for  temperatures  T  >  210°K,  where  £    is  very  nearly  tempera- 

—  wc 

ture  independent.   The  quantity  in  the  first  brace  of  Equa- 
tion (4-2)  indicates  the  water  vapor  emissivity  alone, 
whereas  that  in  the  second  brace  is  the  added  contribution 
by  the  part  of  the  CO   emissivity  transmitted  through  the 
water  vapor  "overlap"  near  15  lim. 
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Equation  (4-2)  was  used  in  the  calculation  of  the  atmos- 
pheric emissivities  over  any  limits  (6,  £i0 )  ,  (4,  10),  ...  , 
(1,  10)  for  which  running  mass  values  U(6,,  10),  C(6,  10), 
etc.  were  available  from  the  radiative  soundings.   Hence 
Equation  (4-2)  is  extremely  versatile  in  its  applicability. 
In  order  to  complete  the  integration  from  level  k  =  1  to 
the  origin  (B  =  0,  £  =  1.0),  the  U(T,  ),  C'(T  )  soundings 
were  taken  as  isothermal  at  T  =  T   over  the  layer  (0,  1). 
However,  it  is  still  necessary,  as  Sasamori  points  out,  to 
determine  the  hatched  areas  on  the  Yamamoto  diagrams  for 
the  temperature-dependent  emissivity  £   |U(0,  10),  T..  ]  . 
This  was  computed  after  Sasamori  as 

I        [0(0,  10),  Tl  =  8.34T  '353  1O*10  U(°'  10)  "  '44 
wc  1  1 

*U(0,  10)"'°3455"'705  *[8.0/(.353  log1Q  U(0,  10)  +  3.56)] 

+  .07262/  (1.  -  .62556[U(0,  10)  +  .02JS6]     ) 

*[log1Q  C(0,  10)  +  1.064]  |  (4-3) 

The  final  term  in  the  braces  of  Equation  (4-3)  involves 
essentially  the  same  form  of  the  CO   correction  to  the 
total  emissivity  as  appeared  in  Equation;  (4-2).   The  first 
term  in  Equation  (4-3)  is,  however,  temperature-dependent, 
as  well  as  dependent  upon  the  summed  absorber  masses  rela- 
tive to  the  reference  level  (taken  here  as  k  =  10) .   It  is 
clear  that  Equation  (4-3)  is  applicable  with  appropriate 
final  parameters  U(0,  k),  C(0,  k),  and  T    to  the 
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determination  of  the  curvilinear  portion  of  the  downward 
flux  contribution  through  level  k  for  all  of  the  reference 
level  cases  to  be  considered  below,  e.g.,  k  =  2,  6. 

B.   NET  FLUX  F   ,  WITH  CLOUDY  SKIES 

The  radiative  soundings  have  been  evaluated  in  the 
form  of  the  parameters  U(k,  10),  C(k,  10),  T  ,  level  for 
level.   In  addition  there  is  listed  for  each  radiative 
sounding  at  every  grid-point  the  cloud-cover  parameters 
CL(1)  and  CL(2),  which  in  general  are  not  both  zero.   The 
grid  area  may  then  be  visualized  as  comprising  the  areal 
fractions  W(0,  0),  W(l,  1),  W(l,  0),  and  W(0,  1)  defined 
by  Equation  (3-12),  of  the  following  cases! 

(a)  totally  clear  skies,  (case  0,  0) 

(b)  overcast  in  both  layers,  (case  1,  1) 

(c)  overcast  in  upper  layer  only,  (case  1,  0) 

(d)  overcast  in  lower  layer  only,  (case  0,  1) 

The  resultant  radiative  effect  of  the  true  cloud  cover 
over  the  grid  area  is  then  reconstructed  by  multiplying 
the  results  of  the  respective  cases  by  the  appropriate 
weight  factors.   It  is  expedient  to  discuss  first  the 

totally  clear  sky  case  (Figure  7(a)). 

* 
1 .   Net  Flux  F    with  Clear  Skies 


If  the  hatched  area  F+    of  Figure  7(a)  is  evalua- 
ted using  the  trapezoidal  summation  rule,  then  the  net 
flux  F    becomes 
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Figure  7(a). 
(case (0,0)  )  . 


Terrestrial  net  flux  F    with  clear  skies, 
The  unhatched  area  depicts 
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F    =  B 
10      10 


.5{e    (8,10)  (B.  -B_)+(B  -B,)  [E    (6,10) 
wc  10   8     8   6    wc 


+  £   (8,10)]  +  [E    (4,10)+£   (6,  10) ]  (B-B . )  +  (B  -B_ )  [£   (2,10) 
wc  wc  wc  6   4     4   2    wc 

+  e   (4,i0)]  +  [e   (i,lo)+e   (2,10)] (b.-b. )+e  •[  (0,10) ,T. ]B. } 

wc  wc  wc  2   1    wc  11 

Bk  =  1.170403  *  10~7  t£  ly(day)"1  (4-4) 

Here  B,   is  the  Stephen-Bol tzman  blackbody  flux  at  T,  . 

* 

2.   Net  Flux  F    with  Overcast  Clouds  in  Both  Layers 

Here  Figure  7 (b)  is  converted  into  a  numerical 
form  using  the  Sasamori  type  emissivities  [Equations  (4-2), 
(4-3)]  in  the  quadrature  process.   The  result  is 


Fio(1'  1]   =  (Bio  "  V[1- 


.5e    (9,  10)] 
wc 


(4-5) 


where  B   is  defined  in  Equation  (4-4). 

* 

3.   Net  Flux  F  with  Overcast  Clouds  in  Upper  Layer 
Only 

Here  Figure  7(c)  is  converted  into  a  numerical 

form  by  the  quadrature  process  after  computing  £    using 

Equations  (4-2)  and  (4-3).   The  result  is 


pio{1'  0)  =  (Bio-V  {1 


5[£    (8,  10) (B.  -B,) 
wc  10   6 


+[e    (8,  10)+C    (6,  10)] (B-B,)] }    (4-6) 
wc  wc  8   6 


The  rationale  of  Equation  (4-6)  is  that  the  downward  flux 
is  essentially  B   plus  an  increment  which  is  in  the  upper 
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Figure  7(b).   Terrestrial  net  flux 
and  low  level  overcast 
unhatched  area  depicts 


F    with  both  high 
and  low  level  overcast  cloud  decks,     (case  (1,1)).   The 
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e<u) 


Figure  7(c).   Terrestrial  net  flux  F    with  high  level 
overcast  clouds  in  the  layer  k=4  to  k=6  only,   (case  (1,0)) 
The  unhatched  area  depicts 
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right  corner  of  Figure  7(c),  given  by  the  summed  combina- 
tion of  the  products  of  the  form 


.5  £    *  A  B 
wc 


as  these  appear  in  the  brackets  in  Equation  (4-6). 

4.   Net  Flux  F    with  Overcast  in  Lower  Layer  Only 

Here  the  numerical  result  follows  immediately  from 
the  graphical  depiction  of  Figure  7(d).  This  case  becomes 
a  replicate  of  the  result  of  Section  IV(B)2,  with 


^o10'11  =  (Bio"  V[1-  -.5ewc<9,lOH 


(4-7) 


5.   Composite  F    Calculations 

Here  we  use  the  principle  enunciated  in  Section  IV 
(B),  and  express  F   [CL(1),  CL(2)]  as  the  properly  cloud- 
weighted  sum  of  the  respective  results  of  cases  (a),  (b), 
(c),  and  (d).   The  result  must  have  the  form  [see  Equations 
(3-13)  and  (3-26)] 


F*Q  =  W(0,0)  F*Q(0,0)  +  W(l,l)  F*0U,1) 

+  W(1,0)  F*   (1,0)  +  W(0,1)  F*   (0,1)         (4-8) 


10 


10 


with  the  cloud-weighted  fractions,  W(0,0),  W(l,l),  W(1,0) 
and  W(0,1)  defined  in  Section  III(B)3.  The  composite  re- 
sult simplifies  to 


F1Q  =   [1-CL(2) ] 


{  (Bio" 


5  Cwc(8'10)   (B10-  B8) 


+  [£wc(8'10)  +  Cwc(6'10)1   (B8  "  V 


}    * 
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Figure  7(d). 

cast  clouds 

The  unhatched  area  depicts 


Terrestrial  net  flux  F    with  low  level  over- 
cast clouds  in  the  layer  k=8  to  k=9  only,  (case  (0,1)). 
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{  V-5  [' 


+<1-CL(2)) (l-CL(l))^  B6-.5  | (£wc(6,10>+ewc(4,10) ) (B6-B4) 
+[£wc(4,l0)+ewc(2,i0)] (B4-B2)+[ewc(2f10)+ewc(lfO)] (b2-b1) 

+  ewcC(0,10)  ,TX)J  I    +CL(2)  I  (B1Q-B9)  [l--.5ewo(9,10)]  } 

(4-9) 

* 
Note  that  Equation  (4-9)  reduces  to  F   (O,0)  when  CL(1)  = 

CL(2)  =  0,  and  to  F   (1,1)  when  CL(1)  =  CL  ( 2 )  =  1.0.   Simi- 

* 
larly  the  other  two  subcases  of  F    may  be  recovered  by 

setting  up  the  weights  corresponding  to  CL  ( 1 )  =  1.0., 
CL(2)  =  0.0  and  CL ( 1 )  =  0.0,  CL(2)  =  1.0„  respectively. 

This  procedure  is  similarly  extended  to  levels  k=6 
and  k=2.   The  results  are  as  listed  in  thve  following  sub- 
sections . 


C.   NET  FLUX  F^  WITH  CLOUDY  SKIES 

D 

Formulation  of  expressions  for  the  unvhatched  areas  of 
Figures  8(a),   (b),   (c)  and  (d),  with  appropriate  weight 
factors  W(0,0),  W(l,l),  W(1,0),  and  W(0,1)  applied,  when 
summed  lead  at  once  to  the  following  result: 


F.  =  (l-CL(l))  ^B  -.5[£    (6,8)  (B  -B,)+£    ( 4  ,  6 )  (B  -B.  ) 
6  (8       wc         8   6    wc         6   4 

+  [E    (4,6)+£    (2,6)  ]  (B  -B.)+[e  r.(2,6)+£    (1 ,  6)  ]  (B.-B.  ) 

wc         wc  4   2     wc         wc  2   1 


+  C    (  (0,6)  ,T 
wc 


1)B1]} 


+  (l-CL(l)  )  (1-CL(2)  ) 


{(Bio- 


B0) [1- 


.5(C    (6,8)+C    (6,10))]!+  CL(1)  /(BQ-BC)  [l-.5(£    (6,8))]1 
wc         wc  )  (86  wc         ) 

+  CL(1)  (1-CL(2)  )   C(B,n-B_)  (l-.5[6    (6,8)+C    (6,10)])} 

I   10   8  wc         wc         ) 

(4-10) 
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(case  (0,0) ) 


Terrestrial  net  flux  F  ,  with  clear  skies 
The  unhatched  area  depicts 
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Figure  8(b).   Terrestrial  net  flux  F  ,  with  both  high  and 
low  level  overcast  layers  (case  (1,1)).   The  unhatched 
area  depicts 
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Figure  8(c).   Terrestrial  net  flux  F 
overcast  clouds  in  the  layer  k=4  to  i 
The  unhatched  area  depicts 


,  with  high  level 
=6  only  (case  (1,0)) 
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Figure  8(d).   Terrestrial  net  flux  F  ,  with  low  level 


overcast  clouds  in  the  layer  k=8  to 
The  unhatched  area  depicts 


■9  only  (case  (0,1)  ) 
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Equation  (4-10)  reduces  to  the  properly  formulated  net 
fluxes  for  cloud  cover  cases  (0,0),  (1,1),  (1,0),  and  (0,1) 
as  these  are  depicted  in  Figures  8(a),(b),(c)  and  (d). 

* 
D.   NET  FLUX  F   WITH  CLOUDY  SKIES 

Similar  to  F, ^  and  F-  the  net  flux  at  level  k=2  can  be 
10       6 

calculated  using  the  trapezoidal  rule.   The  net  flux  F 
is  schematically  depicted  by  the  unhatched  areas  of  Figures 
9(a) , (b) , (c) ,  and  (d) .   After  these  areas  have  been  calcu- 
lated and  weighted  properly  as  described  in  the  preceding 
two  sections  the  following  formula  results 

F*=(1-CL(1) ) {B  -.5 [e    (2,4) (B  -B.)+(e    (2,4)+e    (2,6)) 
2  8      wc         4   2     wc        wc 

(B,-B.)  +  [E    (2,6)+£    (2,8)  ]  (B  -B,)+E    (1,2)(B  -B.) 
6   4     wc        wc  8   6    wc         2   1 

+  £    (  (0,2)  ,Tn  )  (B.  )  ]  }  +(1-CL(1)  )  (1-CL(2)  ) 
wc  1     1 

{(B.  -BQ) [l-.5(£   (2,8)+e   (2,10)]} 
10   8  wc        wc 

+  CL(1)   {B.  -.5[E    (2,4) (B  -B0)+£    (1,2)(B  -B_) 
4        wc         4   2    wc         2   1 

+  £    (  (0,2)  ,Tn  )B.  ]  }  (4-11) 

wc  11 

This  equation  also  reduces  to  the  net  fluxes  for  the  four 
cloud  cover  cases  (0,0),  (1,1),   (1,0),  and  (0,1)  when  the 
appropriate  weight  factors  for  the  simplified  cloud  covers 
are  set  into  Equation  (4-11). 
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Figure  9(a).   Terrestrial  net  flux  F  ,  with  clear  skies 
(case  (0,0)).   The  unhatched  area  depicts 
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CO)) 


Figure  9(b)0   Terrestrial  net  flux  F  ,  with  both  high  and 
low  level  overcast  cloud  decks  (case  (1,1)).   The  unhatched 
area  depicts 
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C(U) 


Figure  9(c).   Terrestrial  net  flux  with  high  level  over- 
cast clouds  in  the  layer  k=4  to  k=6  only  (case  (1,0)). 
The  unhatched  area  depicts 
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Figure  9(d)  „   Terrestrial  net  flux  F  ,  with  lov;  level 
overcast  clouds  in  the  lower  layer,  K=8  to  k  =  9,  only 
(case  (0,1)).   The  unhatched  area  depicts 
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E.   APPLICATION  TO  HEAT  BALANCE  COMPUTATIONS 

In  computing  the  heat  balance  of  the  atmosphere,  the 
layers  (2,6)  and  (6,10)  were  considered  to  be  undergoing 
solar  heating  (Section  III).   Now  we  require  the  long- 
wave cooling  effects  both  instantaneously  and  averaged 
over  the  day  (Section  V).   The  long-wave  cooling  in  (2,6) 
is  symbolized  by  FU26  which  is  defined  by 


*     * 
FU26  =  F„  -  F 
2      6 


Likewise  that  in  (6,10)  is  symbolized  by  FU610  where 


*     * 

FU610  =  F   -  F,_ 
6      10 


is  the  flux-divergence  in  the  lower  half  of  the  atmosphere. 

For  computations  of  the  radiative  "balance"  at  the  tropo- 

* 
pause  only  F   of  Equation  (4-11)  is  required  along  with 

the  net  usable  insolation.   At  the  surface,  a  radiative 

balance  is  computed  at  each  gridpoint  involving  the  absorbed 

* 
insolation  and  the  net  flux  F    of  Equation  (4-9),  the 

latter  considered  a  negative  contribution.   As  will  be 
seen  in  Section  V,  the  long  wave  flux  divergences  FU26  and 
FU610  both  act  as  destabilizing  influences  in  their  respec- 
tive layers . 

F.    STATISTICAL  RESULTS 

1.   Downv/ard  Flux  with  Clear  Skies 


Brunt  (1932)  demonstrated  that  the  nighttime  radia- 
tive measurements  at  the  earth  were  expressible  in  a 
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statistical  form 


Fd  =  Bio(a  +  b/*  > 


(4-12) 


where  B    is  the  Stephen-Boltzman  black  body  radiation 

4 
(B    =  St*T    )  at  the  surface,  e  is  the  water-vapor 

pressure  in  mb ,  and  a  and  b  are  the  regression  coeffi- 


cients, and  F,  =  B,.  -  F,,..   The  regression  coefficients 
d     10     10 

referred  to  by  Brunt  were  shown  to  be  subject  to  wide 
variations  depending  on  B    and  /e. 

In  this  study  a  statistical  test  of  Equation  (4-12) 
was  carried  out  by  linear  regression  using  a  program  from 
the  BMD  set  of  statistical  programs  (Dixon,  1973) o   The 


predictand  selected  for  these  computations  was  F  ,  which 

included  the  downward  contribution  of  CO   as  well  as  that 

of  water  vapor  on  the  calculation  of  the  net  flux  at  level 
10.   The  statistical  equation  that  resulted  was 


Fio  =  B10(-30183  " 


01314  /e) ,  R  =  .995 


(4-13) 


The  correlation  coefficient  of  Equation  (4-13)  shows  that 

B    and  B    /e  are  good  predictors  of  F    and  therefore 

of  F  .   Using  the  values  of  a  and  b  in  Equation  (4-13)  and 

solving  for  F,  we  obtain 
a 


or 


Fd  =  Bio 


(1-a)  -  b/i 


F0  =  B,rt(. 69817  +  .01314  /e) 
a     10 


(4-14) 


69 


The  high  value  of  the  correlation  coefficient  indi- 
cates that  Equation  (4-13)  may  be  regarded  as  an  oceanic 

version  of  the  Brunt  F  -equation.   Note  however  that  most 

d 

of  the  dependence  of  F   is  with  B    and  that  the  variation 

*  /— 

of  F    with  B   /e  is  relatively  slight.   This  is  as  to  be 

expected  over  the  ocean,  where  nearly  constant  relative 
humidities  might  well  exist  over  large  areas.   The  mean 
values  of  F   and  B    for  clear-sky  cases  are  as  listed  be- 
low 


F,  =  637.09  ly  day"1 
d 

B   =  843.47  ly  day"1 

so  that  F,  averages  75.5  percent  of  Biri.   These  statistical 
d  10 

averages  are  based  upon  the  57  grid-point  soundings  in  the 

Northern  Hemisphere  only  (those  in  the  Southern  Hemisphere 

have  been  deleted  in  order  to  eliminate  tropical  bias) . 

2 .   Statistical  Relationship  between  Cloudy  Sky  and 
Clear  Sky  Cases 

Based  upon  the  radiative  sounding  parameters  listed 

in  Table  1(b)  and  Equation  (4-9)  it  was  possible  to  compute 

*  * 

F   [CL (1 ) , CL (2 ) ]  and  F   (0,0)  based  upon  otherwise  identi- 
cal parameters. 

For  the  purpose  of  this  statistical  study  the  ratio 

*  * 

F   [CL  ( 1 )  , CL  (2 )  ] /F   (0 , 0 )  was  formed  as  the  predictand, 

while  an  approximation  to  the  composite  total  opaque  cloud 
cover  at  the  grid-point,  defined  in  Equation  (3-43)  as 

CL  =  CL(1)  +  CL  (2)  (l-CL(l)  ) 
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for  each  cloudy  sky  case  was  used  as  the  predictor.   The 
simple  regression  formula  based  upon  a  best-fit  to  the  57 
Northern  Hemisphere  data  soundings  became 


fT„ (CL)  =  F*  (0,0) (1  - 


10 


10 


7575CL) 


(4-15) 


with  a  regression  coefficient,  R  =  .983.   A  regression 

2 
analysis  using  CL   as  the  predictor  resulted  in  only  a 

very  small  improvement  of  R.   The  hemispheric  mean  values 

of  the  various  parameters  in  Equation  (4-15)  were  as 

follows 


F*  (CL)  =  112.48  ly  day"1 


F*  (0,0)=  206.38  ly  day  1 


10 


(4-16) 


CL  =  0.6103 


Here  the  superior  bar  means  the  sample  mean  was  taken 
over  the  57  Northern  Hemisphere  grid-points. 

The  mean  results  of  (4-16)  indicate  that  the  sur- 
face net  flux  is  decreased  by  the  ratio  (93.9/206.38)  =  .455 
with  CL  =  0.6103.  This  is  to  say  that  the  downward  flux 

F   ._  is  on  the  average  45.5  percent  larger  than  with  the 
d ,  10 

clear-sky  case  for  the  cloud-parameterization  of  Equation 
(2-13)  . 
3 


*  * 

F    (CL)  Relative  to  F   (0,0) 


Here  the  statistical  parameter  sought  is  the  pre- 


dictand 


F2(CL) 
F*  (0,0) 


(4-17) 
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where  both  numerator  and  denominator  are  computed  directly 
from  Equation  (4-11)  with  (a)  the  known  CL(1) ,  CL ( 2 )  used 
in  the  numerator,  and  with  (b)  CL(1)  =  CL(2)  =  0  in  the 
denominator.   The  regression  result  was  cast  in  the  desired 
form 


Y  =  1  -  A1 (CL)  -  A2 (CL) 


with  A  ,  A   to  be  determined.   The  resulting  regression 


equation  was 


Y  =  1  -  ,0742(CL)  -  ,1450(CL) 


(4-18) 


with  a  multiple  correlation  coefficient  R  =  0.9431,  and  a 
standard  error  of  estimate  .0529.   The  mean  values  of 
F2 (CL) ,  F2 (0,  0) ,  were 


F*(CL)  =  453.19  ly  day  1 


F*Q(0,  0)  =  517.70  ly  day"1 
CL  =  0.6803 

The  mean  CL  here  was  based  upon  the  75  sounding-sample 
which  contained  a  larger  proportion  of  tropical  cases  than 
did  those  of  the  preceding  section. 

It  is  to  be  noted  that  the  infrared  loss  to  space 
has  not  been  greatly  impeded  by  the  assumption  of  normal 
cloud  cover  since  no  thick  cirrus  clouds  have  been  in- 
cluded in  the  model.   With  deep  cirrus  the  "greenhouse 
effect"  would  be  increased  for  gridpoint  columns  so  affected 
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For  comparison  with  satellite  long  wave  flux  to 

.  -1 

space,  which  is  usually  given  in  ly  min   ,  the  foregoing 

*  -1 

F   values  should  be  divided  by  1440  min  day 

Comparison  with  satellite  measures  of  long  wave 

* 
flux  should  be  made  with  F  +  rather  than  with  F  ,  which 

includes  a  downward  flux  contribution  from  the  strato- 
sphere.  For  this  purpose,  computations  were  made  with 
F_  of  Equation  (4-11)  but  omitting  terms  involving  both 
levels  "2"  to  "1"  and  "2"  to  "0",  which  are  always  asso- 
ciated with  downward  flux  (see  Figure  9).   The  comparison 
of  the  grid-point  sounding  mean  F  [CL ( 1 )  ,CL  ( 2 ) ]  was  then 
made  with  the  NIMBUS  III  global  average  long  wave  flux 
to  space  for  the  appropriate  date  (after  Raschke  et  al, 
1973).   The  comparison  is  as  follows 

F2[CL(1),  CL(2)]=  .407  ly  min"1 
long-wave  flux  to  space=  .352  ly  min     (after  Raschke) 

The  comparison  is  quite  good  considering  that  the  results 
of  Raschke  et  al  are  Northern  Hemisphere  means  for  a  two- 
week  period  in  April  1969,  whereas  the  results  of  this 
study,  .407  ly  min   ,  is  essentially  based  upon  a  mid- 
Pacific  cros s- sec tion  on  April  25th,  which  tends  to  be 
warmer  than  the  globe  as  a  whole. 
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V.   MERIDIONAL  CROSS-SECTIONAL  DEPICTION  OF  RESULTS 

In  the  computer  programming  of  this  model  many  outputs 
have  been  generated,  so,  for  ease  of  examination  of  the 
output,  an  atmospheric  cross-section  along  each  of  the 
three  meridians  studied  was  developed.   This  form  of  pre- 
sentation was  designed  to  give  both  the  incoming  and  the 
outgoing  radiations  at  the  selected  data  levels.   In  addi- 
tion, radiation  balances  were  computed  within  key  atmos- 
pheric layers  and  at  the  earth's  surface  at  each  gridpoint. 
In  order  to  indicate  the  relevance  of  the  cloud-amounts 
CL(1),CL(2)  and  CL  to  the  resulting  balances,  these  values 
were  also  plotted  at  appropriate  levels  (Figure  11,  in 
parentheses)  at  gridpoints.   These  cross-sections  were 
then  interpolated  to  give  the  above  results  at  five  degree 
latitude  gridpoints  along  the  selected  meridians.   This 
was  done  for  ease  in  calculation  of  a  composite  value  of 
each  of  the  computed  radiative  quantities  at  identical 
latitudes  between  20°S  and  65°N,  wherever  possible. 
Lines  1  and  3  extended  from  20°S  to  55°N  and  had  identical 
latitude  grid-points.   Line  2  extended  from  the  equator  to 
65°N  and  therefore  permitted  extension  of  the  cross-section 
to  65°N.   At  latitudes  south  of  the  equator  only  two  values 
were  available  for  the  composite  computations  and  at  lati- 
tudes higher  than  55°N,  only  one  value  was  available. 
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Between  these  extremes  all  three  meridional  cross-sections 
were  used  in  the  calculation  of  the  composite  case.   Fig- 
ure 11(a)  and  (b)  depicts  the  composite  cross-section  with 
Figure  10  providing  the  key  for  locating  the  desired  values. 

In  the  cross-section  to  be  depicted  the  insolation 
rates  at  level  k=2,  (QAVG)  are  given  as  averages  over  the 
24-hour  day  (25  April  1973).   This  was  done  by  employing 


Equation  (3-1),  but  replacing  Cos  Z  by  Cos  Z  where 


Cos  Z  =  [h  Sin<f>  Sin6  +  Costf)  Cos6  Sin  h]/TT       (5-1) 

and    h  =  hour  angle  =  Arc  Cos  (-Tarufi  TanS) 
<f>  =  latitude 
6  =  declination  (Smithsonian  Tables  (List,  1963)) 

Here  the  double  bar  over  Cos  Z  represents  a  time  mean  over 
a  24  hour  period,  that  is  for  a  time  period  compatible 
with  the  long-wave  flux  transfers. 

At  each  latitude  a  value  of  the  instantaneous  absorp- 
tivities,  reflectivities  or  transmissivi ties  was  known  for 
the  time  0000GMT,  and  for  the  purpose  of  the  cross-section 
presented  here,  these  were  assumed  constant  even  though 
the  insolation  rates  changed  over  the  day.   The  IR  trans- 
fer rates  and  their  flux  divergences  were  assumed  to  be 
constant  both  for  the  day  and  the  composite  grid-point  as 
determined  from  their  0000GMT  soundings.   Thus  it  was 
possible  to  compute: 
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a)  the  earth-atmosphere  radiative  balance  at  level 
k=2  by  the  formula 

* 
QTROP  =  QAVG  -  (REFF  +  F  )  (5-2) 

b)  the  radiative  balance  in  the  layer  (2,  6)  as  given 
by 

NET26  =  ABS26  +  FU26  (5-3) 

c)  the  radiative  balance  in  the  layer  (6,  10)  as  given 
by 

NET610  =  ABS610  +  FU610  (5-4) 

d)  the  radiative  balance  at  the  earth's  surface  as 
given  by 

NETSFC  =  ABSFC  +  F  (5-5) 

In  addition  to  these  mean  grid-point  items  the  cloud  covers 
CL(1),  CL(2)  and  (CL)  in  parentheses  are  depicted  in  the 
applicable  layer  after  averaging  to  the  grid-point.   At  the 
tropopause,  and  at  the  earth's  surface  a  global  average  of 
radiative  balance  or  net  was  obtained  by  cosine  weighting 
as  follows 


]£  Q  (<!>!>  Cos  <J> 
Q  =  —  (5-6) 


Cos  d>  . 
i 


where  Q[<p.)   is  the  radiative  net  value  under  consideration 

at  latitude  it.   .   Some  of  these  values  are  contrasted  with 
l 
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the  other  climatological  results  cited  below 


1)  global  balance  at 
the  tropopause 

2)  global  balance  at 
the  surface 


This  Study    Other  Data 


048 


166 


061,  Raschke 


.140,  Budyko 


The  outgoing  long  wave  radiation  in  these  computations  was 
also  in  excellent  agreement  with  that  of  Raschke  et  al 
(1973)  as  noted  in  Section  IV.   The  global  balance  of  total 
radiation  at  the  tropopause  is  less  than  that  reported  by 
Raschke  et  al  for  the  same  general  time  of  year.   However, 
it  was  noted  that  this  result  primarily  reflects  the 
smaller  value  of  planetary  albedo  (ALB  =  .300)  reported 
from  satellite  observations. 

Another  set  of  results  for  comparison  with  those  calcu- 
lated here  are  due  to  Budyko  (1956).   Budyko  calculated 
the  radiation  balance  at  the  surface  for  ten-degree  lati- 
tude bands.   In  his  study,  these  values  were  taken  for  the 
balance  at  mid-points  of  these  bands  and  compared  with  the 
computed  output  from  the  present  model.   The  following 
table  shows  this  comparison. 


5°N 
15°N 
25°N 
35°N 
45°N 
55°N 


This  Study 
.  158 
.  205 
.  236 
.  201 
.119 
.144 


Budyko 

(Longitude  180°) 
.  203 
.  209 
.  190 
.  148 
.  103 
.064 
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-2      -1 
All  of  the  above  values  are  in  cal  cm    min   .   These 

comparisons  would  have  been  in  closer  agreement  if  the 

results  of  this  study  had  been  in  the  form  of  an  annual 

average  as  Budyko's  results  were. 
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VI .   CONCLUSIONS 

Using  three  meridional  lines  of  soundings  over  the 
Pacific  Ocean  a  radiative  model  for  inclusion  in  the 
FNWC  numerical  weather  prediction  scheme  was  evaluated. 
These  three  meridians  of  data-soundings  were  used  to 
compute  a  space-time  average  of  the  radia-ttional  para- 
meters for  one  day  over  the  Pacific  Ocean.   This  space- 
time  average  agreed  fairly  well  with  values  derived  from 
other  studies,  including  satellite  measurements.   There 
were  some  minor  differences  which  could  be  ascribed  to 
the  fact  that  the  data  used  in  this  experiment  was  a  one 
day  space-time  average  while  the  comparison  data  was  an 
average  of  several  weeks  for  the  satellite  data,  and 
over  years  in  the  case  of  surface  climatological  data. 

Differences  could  also  be  induced  by  the  parameteri- 
zation of  cloud  amounts  as  well  as  by  the  treatment  of 
cloud  reflectivity  and  absorptivity  for  solar  wavelengths 
In  this  model,  only  two  layers  of  clouds  were  considered 
and  the  same  empirical  values  of  the  reflectivities  and 
absorptivi t ies  were  used  for  all  latitudes,  which  could 
cause  discrepancies  when  compared  with  satellite  data. 

Although  the  results  of  this  data  treatment  agree 
fairly  well  with  climatological  comparisons  cited  at  the 
end  of  Section  V,  refinements  may  be  possible  in  the 
present  model.   This  model  should  be  tested  using  data 
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at  different  seasons  of  the  year  and  over  varying 
locales.   With  a  full  year  of  sampling,  refinements  in 
the  cloud  parameterization,  atmospheric  absorption 
values  and  reflection  treatments  could  be  developed 
for  accurate  use  in  the  FNWC  predictive  model. 
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